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' Abstract 

a \ 

We utilise the fact that the CCFM equation in the single loop approximation can be 
diagonalised by the Fourier-Bessel transform. The analytic solution of the CCFM equa- 
tion for the moments f u (b, Q) of the scale dependent gluon distribution is obtained, where 
b is the transverse coordinate conjugate to the transverse momentum of the gluon. The 
unintegrated gluon distributions obtained from this solution are analysed. It is shown 
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unintegrated gluon distributions in terms of the integrated ones. The corresponding ap- 



proximate expressions for the unintegrated gluon distribution are compared with exact 



solution of the CCFM equation in the single loop approximation. 



1 Introduction 



The basic, universal quantities of the QCD improved parton model are the scale dependent 
parton distributions, like the gluon distribution g(x, Q 2 ), where x denotes the momentum frac- 
tion. The (integrated) parton distributions can be related to the less inclusive distributions 
f(x, Qt, Q) unintegrated over transverse momentum Q t of the parton. Those uninegrated dis- 
tributions are often needed in less inclusive measurements which are sensitive to the transverse 
momentum of the parton [|I|] - ||. 

The unintegrated, scale dependent distributions are described in QCD by the Catani- 
Ciafaloni-Fiorani-Marchesini (CCFM) equation - |17| based upon quantum coherence which 



implies angular ordering |fL8|j . Its very nice feature is the fact that it embodies both the DGLAP 
and BFKL evolutions at low x. In the region of large and moderately small values of x, where the 
small x effects can be neglected the CCFM equation becomes equivalent to the (LO) DGLAP 



evolution. This approximation corresponds to the so called 'single loop' approximation 0, 10 



The CCFM equation interlocks in a rather complicated way the two relevant scales i.e. the 
transverse momentum Q t of the parton and the hard scale Q. The main purpose of this paper 
is to explore the fact that in the 'single loop' approximation the CCFM equation can be solved 
exactly in the transverse coordinate representation conjugate to the transverse momentum of 
the parton. The unintegrated distributions can then be obtained from the Fourier-Bessel trans- 
form of this solution. Although the 'single loop' approximation neglects small x effects and so it 
is not valid at very small x it can be a reasonable approximation at large and moderately small 
values of x ( x > 0.01 or so), which is certainly the region of phenomenological interest ||, |5|. 
Analytic insight into the exact solution of the CCFM equation in the single loop approximation 
will also make it possible to critically examine and justify approximate formulas linking the 
unintegrated distributions to the integrated ones in the region where the LO DGLAP dynamics 
should be adequate jl], |], 



The content of our paper is as follows: In the next Section we recall the CCFM equation 
for the unintegrated gluon distribution. In Section 3 we discuss the 'single loop' approximation 
of this equation in the transverse coordinate representation. We show that it can be solved 
exactly for the moment function Q), where b denotes the transverse coordinate conjugate 
to the transverse momentu Q t of the gluon. We do also show how the approximate forms 
of this solution expressing the unintegrated distributions in terms of the integrated ones (l|. 
originate from the exact solution. In Section 4 we present numerical results for the unintegrated 
gluon distributions based on the solution of the CCFM equation in the transverse coordinate 
representation. We also confront exact solution with its approximate forms. Finally in Section 
5 we summarise our main results and give our conclusions. 
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2 The CCFM equation 



Parton cascade with angular ordering generates the Catani, Ciafaloni, Fiorani, Marchesini 
(CCFM) equation for the unintegrated, scale dependent gluon distribution f(x,Q t ,Q) in 
the proton, where x, Qt and Q denote the longitudinal momentum fraction carried by the gluon, 
transverse momentum of the gluon and the hard scale respectively. The latter is specified by 
the maximal value of the emission angle. The CCFM equation has the following form: 



/(*, Q t , Q) = f°(x, Q t , Q)+ / ^| f -Q(Q - qz)Q{q - q )^A s (Q, q, 



2N c A NS (Q t , q, z) + + zP gg (z) 



1-z) 



X 



f[-\Q t +(l-z)q\,q) 



(1) 



where A$(Q, q, z) and A^s{Qt, Q, z ) are the Sudakov and non-Sudakov form factors. They are 
given by the following expressions: 



&s(Q,Q,z) = exp 



dp 2 a s fi-Qo/p 



(qz) 2 p 2 2ir Jo 



dZZPgg(z) 



&Ns{Qt,q,z) = exp 



1 dz' [Qt dp 2 2N c a s 



lz Z' J(qz') 2 p 2 2TT 

For simplicity we neglect possible quark contributions. The function P gg (z) is: 

P gg (z) = 2N c [-2 + z(l - z)} 
and corresponds to the non-singular part of the g — > gg splitting function P gg (z) 



P gg (z) = 2N C 



1 1 

+ 



z 1 — z 



+ P gg (z) 



The argument of a s will be specified later. 



(2) 
(3) 

(4) 
(5) 



In principle the CCFM equation has been obtained using only the singular parts of the 
splitting function proportional to 1/z and 1/(1 — z). We add the non singular part P gg (z) to 
the kernel of this equation in order to obtain the complete DGLAP evolution in the 'single 
loop' approximation. The two-dimensional vector q in equation (P is related to the transverse 
momentum q t of the emitted gluon 



Qt = (1~z)Q 



(6) 



The constraint Q > qz reflects the angular ordering and the inhomogeneous term f°(x,Q t ,Q) 
is related to the input non-perturbative gluon distribution. It also contains effects of both the 
Sudakov and non-Sudakov form-factors. 

It should be observed that if the cut-off qz' in the definition of the non-Sudakov form-factor 
is replaced by the fixed cut-off go then the non-Sudakov form-factor reduces to the form-factor 
reflecting the reggeisation of the gluon, i.e.: 

qz' — g -> &ns = exp[(2a G (Q 2 ) - 2)ln(l/z)] 

t dp 2 N c a s 



it P 2 2tt 



The unintegrated scale dependent gluon distribution f(x, Q t , Q) is related in the following 
'standard' way to the conventional (scale dependent) integrated gluon distribution xg(x,Q 2 ): 



xg(x,Q 2 ) = f dQ 2 f(x,Q t ,Q) (7) 

In the 'single loop' approximation of the CCFM equation ([]]) the angular ordering constraint 
Q(Q — qz) is replaced by Q(Q — q) and the non-Sudakov form-factor A NS is set equal to unity 
II |TU|. Equation then reads: 



f(x,Qt,Q)=f°(x,Qt,Q)+[ ~\ I -e(Q-q)e(q-q )^A s (Q,q,z=l)* 

/(f,lQt+(l-«)9l,?)) (8) 

It is useful to 'unfold' the Sudakov form factor in equation @ in order to treat the real emission 
and virtual corrections terms on equal footing. Unfolded CCFM equation in the single loop 
approximation takes the following form: 

f(x, Q t , Q) = f°(x, Q t ) + I — ^9(g - g )^ /' -zP n (z)* 

J 7cq z 2tt Jo z 

9(Q - q)Q(z - x)f \Q t + (1 - z)q\, g)) - zB(Q - q)f(x, Q t , q) 

The inhomogeneous term f°(x, Qt) is equal to the input non-perturbative gluon distribution in 
x and Q t . 



™ c + ^^ ) +zP gg (z) 



(9) 
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3 Transverse coordinate representation of the CCFM 
equation in the single loop approximation 

It can easily be observed that the CCFM equation in the single loop approximation @ can be 
diagonalised by the Fourier-Bessel transform: 

POO _ 

f(x,Q t ,Q)= dbbJ (Q t b)f(x,b,Q) (10) 
J o 

with the function f(x, b, Q) given by: 

f(x,b } Q) = dQ t QMQ t b)f(x,Q t ,Q) (11) 
Jo 

where Jq{u) is the Bessel function. From equations (0) and (pTTf ) we get: 

f(x,b = 0,Q) = ^xg(x,Q 2 ) (12) 

The corresponding equation for f(x,b,Q), which follows from equation (|^) after taking the 
Fourier-Bessel transform of both sides of this equation reads: 

/>,6, Q) = + f d 4^1 f * xPM . 

Jql (7 Z7T JO Z 



^e(z-x)J [bq(l-z)q]f(^,b,q)^-zf{x,b,q)^ (13) 

where we put q 2 as the argument of a s . This choice of scale gives standard LO DGLAP equa- 
tion for the integrated gluon distribution for xg(x, Q 2 ) = 2f(x, 6 = 0, Q). 

In order to solve equation C p~3|) it is useful to introduce the moment function f^(b, Q) 

L(b, Q) = C dxx^fix, b, Q) (14) 
Jo 

Equation ( |T3"D implies the following equation for the moment function f w (b,Q): 



Ub, Q ) = ;» + r %^fl [dzz Pgg {zy 

Jq^ q ZTT JO 

{^- x Jo[6g(l - z)q]f u (b,q) ~ f u Q>,q)} (15) 
The solution of this equation reads: 

M,Q)=fi(b)exp[S u {b,Q)] (16) 



4 



where 



rQ 2 do 2 a (a 2 ) r 1 

S u (b,Q) = I Jr^T / dzzP gg {z){z^J Q [{l - z)bq] - 1} (17) 

Jql q z Z7T JO 

At small values of b (i.e. b « 1/qo) we can neglect b dependence in and set 

£(&) * f"(b = 0) 

We can identify f®(b = 0) with the moment of the input (non-perturbative) integrated 
distribution, i.e. 

f> = 0) = \& (18) 

where 

gl = J dQ 2 t £ dxx^fix, Q t ) (19) 

We note that at b = solution ([16]) reduces to the solution of the DGLAP equation for the 
moment function g^iQ 2 ) of the integrated gluon distribution g(x,Q 2 ), i.e. 

9^Q 2 )= t dxx"g(x,Q 2 ) (20) 
Jo 

To be precise we get 

f(b = 0,Q) = ±gUQ 2 ) 
UQ 2 ) = glexp \ f f dzzP gg (z)[z^ - 1]} (21) 

yJq q ^7T JO J 



It is useful to rearrange solution ([Tq ) as below: 

U(b,Q)=U(b,Q)T g (b,Q) (22) 

where 



Q 2 dq^ajrf} i' 1 
g 2 27r Jo 



/«(&.<?) = OTexp^ / ^T^t 1 L dzzP gg {z)J Q [{l-z)bq}\z^ - 1] \ (23) 



and the Sudakov-like form-factor T g (b,q) is given by: 

T g (b, Q) = exp { jff jT 1 ^P M (,)[J [(1 - z)bq] - 1]} (24) 

In order to obtain more insight into the structure of the unintegrated distribution which 
follows from the CCFM equation in the single loop approximation it is useful to adopt the 
following approximation of the Bessel function: 

J (u)~Q{l-u) (25) 
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Using this approximation in equation (11) we get: 



df(x,b=l/Q t ,Q) 
dQ\ 



f(x,Q t ,Q)^2 " J ^":^^ (26) 



It may be useful to analyse solution ([TB]) using approximation (|^) which gives 

Lib, Q) ~ ^T g (l/b, Q)exp[S r UJ (b 1 Q) + AS^b, Q)] (27) 

where 

r min{i/b 2 ,Q 2 ) da 2 a( a 2 ) r 1 
SL(b,Q) * f / dzzP gg {z)[z^ - 1] (28) 

Jq 2 q £ Z7T JO 



and 



AS:(b,Q) = f Q [ f 1 dzzP^z)^- 1 - 1] (29) 

Ji/b 2 q 2 2n Ji-i (bq) 



h/b 2 q 2 2tt Ji-i/(bq) 
T g (b,Q)~exp(S°(b,Q)) (30) 



where S°(b, Q)) is given by: 



It may be seen that the form-factor T g (l/b,Q) given by equations (^,^) has the structure of 
the Sudakov form-factor. It can also be seen that the factor g®exp[S^(b, Q)} in equation (|2~7f) 
with S^(b, Q) defined by equation ( |28"D can be identified with the moment function g^n 2 ) of the 
integrated gluon distribution at the scale /i 2 = min(l/b 2 ,Q 2 ). Neglecting the term AS^(b,Q) 
we get: 



that gives: 



Mb, Q) ~ T g (b, Q)g u) (mm(l/b i , Q 2 )) (32) 



MQt, Q) - 2 ^ 

d[T 9 (b=l/Q t ,Q) 9uJ (Q 2 )} 



dQ 



(33) 



for Qt < Q, and 



MQt,Q) = o 
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for Q t > Q. Equation fl33|) gives: 



f(x, Q t , Q) - ^-vvw^a^j (34) 



d[T g (b=l/Q u Q)xg(x,Q 2 )} 
dQl 

It should be noted that equation ( |33"D correctly reproduces the double logarithmic effects in the 
region Q t « Q 0. The formalism presented above is similar to that used for the description 
of the pt distributions in (for instance) Drell-Yan process (see e.g. [|Hfl). 

Taking approximately into account the remaining contribution in equation (^) gives: 



Derivation of equation ( p5]) is sketched in the Appendix. 



After replacement Q 2 /(l — z) 2 — > Q 2 in the argument of a s and in the gluon distribution 
g(x/z, /i 2 ) which introduces subleading effects, expression (|35| ) coincides with the representation 
used in ref. @ (modulo subleading terms in the definition of the Sudakov form-factor): 



T 9 (b=l/Q t ,Q) fi~~Qt/Q a s (Q 2 



f(x, Qt, Q) * gv n " / dzP gg (z)^^e(z -x)-g{ -, Qi (36) 

LJf Jo An z \z J 

4 Numerical results 

In the previous Section we have shown that the CCFM equation in the single loop approxima- 
tion can be solved analytically in the h space, where b is the transverse coordinate conjugate to 
the transverse momentum Q t of the gluon. We have also indicated approximations which make 
it possible to relate the unintegrated distributions to the unintegrated ones. In this Section we 
present results of the numerical analysis of the exact solution of the CCFM equation utilising 
its diagonalisation in the transverse coordinate representation. We shall also confront this exact 
solution with the approximate expressions defined by equations fl34|) and (j36|). 

To this aim we solved equation (|13|) for the distribution f(x, b, Q) and computed the unin- 
tegrated distribution f(x,Q t ,Q) from the Fourier-Bessel transform: 

POO _ 

f(x,Q t ,Q)= dbbJ (bQ t )f(x,b,Q) (37) 
Jo 

We started from the input distribution f°(x,b) 



f°(x,b) = 9 -^exp(-b 2 q 2 Q /A) 
7 



g (x) = 3(l-x) 5 (38) 

where we have set go = lGeV. In Fig. 1 we plot the function Q 2 f(x, Q t , Q) as the function of 
Qt for Q 2 = lOOGeV 2 and for two valus of x, x = 0.01 (upper curve) and x = 0.1 (lower curve). 
In Figure 2 and 3 we compare those exact solution with approximate expressions ([34]) and (|36D . 
We find that equation fl36|) gives somehow better approximation of the exact solution except for 
the 'end points' Q\ ~ Q\ and Q\ ~ Q 2 . The simple formula ( [34] ) is a reasonable approximation 
of the exact solution for small values of Qt- It may however give negative contribution at large 
x(x ~ 0.1) and large Q\ . 



Summary and conclusions 

In this paper we have utilised the transverse coordinate representation of the CCFM equation 
in order to get an analytical insight into its solution. The transverse coordinate representation 
has been widely used for the discussion of the soft gluon resummation effects in the transverse 



momentum distribution of Drell - Yan pair etc., ||19|| . In our paper we have utilised the fact that 



this representation diagonalises the CCFM equation in the single loop approximation and can 
be very helpful for obtaining unintegrated parton distribution satisfying the CCFM equation in 
this approximation. We have shown that the CCFM equatiion in the single loop approximation 
can be solved analytically for the moment function Q), where b is the transverse coordinate 
conjugate to the transverse momentum of the gluon. We have also confronted the unintegrated 
gluon distributions with approximate expressions which were discussed in the literature. The 
single loop approximation neglects small x effects in the CCFM equation and, in particular, 
it neglects virtual corrections responsible for the non-Sudakov form-factor. This form-factor 
generates contributions which are no longer diagonal in the b space and so the merit of using 
this representation beyond the single-loop approximation is less apparent. However, in the 
leading ln(l/x) approximation at small x the CCFM equation reduces to the BFKL equation 
with no scale dependence and the kernel of the BFKL equation in the b space is the same as 
the BFKL kernel in the (transverse) momentum space. One can expect that the transverse 
coordinate representation of the CCFM equation may eventually appear to be helpful beyond 
the single loop and BFKL approximations. 
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Appendix 

In this Appendix we derive equation fl35|) . To this aim we start from the following improved 
approximation of the derivative df w (b = l/Q t ,Q)/dQf: 

2 df u (b = l/QuQ) „ d[T g (b = l/Q t ,Q)g u (Q*)] + = , 0AS*(6 = 1/Q t> Q) 



dQi 



dQi 



(39) 



where the function AS r (b = l/Q t ,Q) is given by equation fl2D|). Using equations (|3Tf ) and 
we get for Q t < Q: 

JfUb = i/Qt,Q) 



T g (b = l/Q t ,Q) 



Q 



l-Qt/Q 



dZZPgg(Z 



a. 



2tt 



9u{Qt) + -JTZTTW^ 



dln(Q 



+ 



T g (b=l/Q t ,Q) 



Q 



2« w ^ ui-j ^ Pgg(z )[^-i _ i] _ / ^^pz 2 p 99 ( 2) [^-i _ l] 



2tt 



1 a.^ 
2tt 



(40) 



Taking into account the DGLAP evolution equation: 



2 dgM^) « s (/i 2 ) f 1 , , u u _! , , 2 N 



we get: 



2tt Jo 



(41) 



T g (6=l/Q f ,Q) 



i-Q*/Q a s (tMv) / 

(fa ^ z) J zP gg {z)z^ 
2tt 



(42) 



Taking the inverse Mellin transform of both sides of equation (JE3) we get equation (|35|). 
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Figure 1: Function Q 2 f(x,Q t ,Q), where f(x,Q t ,Q) is the unintegrated gluon distribution 
obtained from the exact solution of the CCFM equation in the single loop approximation, 
plotted as the function of the transverse momentum Q t of the gluon for Q 2 = lOOGeV 2 . The 
upper and lower curves correspond to x — 0.01 and x — 0.1 respectively. The transverse 
momentum Q t is in GeV. 
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Figure 2: Function Q 2 f(x,Q t ,Q), where f(x,Q t ,Q) is the unintegrated gluon distribution 
obtained from the exact solution of the CCFM equation in the single loop approximation, 
plotted as the function of the transverse momentum Q t of the gluon for Q 2 = lOOGeV 2 and 
x = 0.01. The solid curve corresponds to f(x,Qt,Q) obtained from exact solution of the 
CCFM equation in the single loop approximation, while the short dashed and long dashed 
curves correspond to approximate expressions for f(x, Q t , Q) given by equations (34) and (36) 
respectively. The transverse momentum Q t is in GeV. 
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Figure 3: Function Q 2 f(x,Q t ,Q), where f(x,Q t ,Q) is the unintegrated gluon distribution 
obtained from the exact solution of the CCFM equation in the single loop approximation, 
plotted as the function of the transverse momentum Q t of the gluon for Q 2 = lOOGeV 2 and 
x = 0.1. The solid curve corresponds to f(x,Qt,Q) obtained from the exact solution of the 
CCFM equation in the single loop approximation, while the short dashed and long dashed 
curves correspond to approximate expressions for f(x, Q t , Q) given by equations (34) and (36) 
respectively. The transverse momentum Q t is in GeV. 
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